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ABSTRACT 

The differential migration of two planets due to planet-disk interaction can result in 
capture into the 2:1 eccentricity- type mean- motion resonances. Both the sequence of 2:1 
eccentricity resonances that the system is driven through by continued migration and 
the possibility of a subsequent capture into the 4:2 inclination resonances are sensitive 
to the migration rate within the range expected for type II migration due to planet-disk 
interaction. If the migration rate is fast, the resonant pa ir can evol ve into a family of 



2:1 eccentricity resonances different from those found by iLed (120041 ) . This new family 
has outer orbital eccentricity ei > 0.4-0.5, asymmetric librations of both eccentricity 
resonance variables, and orbits that intersect if they are exactly coplanar. Although 
this family exists for an inner-to-outer planet mass ratio m^/m o > 0.2, it is possible 



to evo lve into this family by fast migration only for m\/m2 > 2. iThommes fc: Lissauer 



(|2003l ) have found that a capture into the 4:2 inclination resonances is possible only 
for mi/m-2 < 2. We show that this capture is also possible for mxjm^ > 2 if the 
migration rate is slightly slower than that adopted by Thommes h, Lissauer. There is 
significant theoretical uncertainty in both the sign and the magnitude of the net effect 
of planet-disk interaction on the orbital eccentricity of a planet. If the eccentricity is 
damped on a timescale comparable to or shorter than the migration timescale, e2 may 
not be able to reach the values needed to enter either the new 2:1 eccentricity resonances 
or the 4:2 inclination resonances. Thus, if future observations of extrasolar planetary 
systems were to reveal certain combinations of mass ratio and resonant configuration, 
they would place a constraint on the strength of eccentricity damping during migration, 
as well as on the rate of the migration itself. 



1. INTRODUCTION 

Extrasolar planet searches have to date yielded about 33 systems with multiple planets, and at 
least 8 of these systems have a pair of planets known or suspected to be in mean-motion resonances. 
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It is well established that the outer two planets in the GJ 876 system are deep in 2:1 resonances, 
with the retrograde peria pse precessions induced by the 2:1 resonances having been observed for 



more than one full period (IMarcv et al 



Lee & Peald 120021 : lLaughlin et al 



2005 



2001 



Rivera et al 



Laug hlin &: Chambersll200ll : iRivera &: Lissauerll2001 



2005J). In the GJ 876 system, both of the 



lowest order, eccentricity- type mean-motion resonance variables 

(9i = Ai - 2A 2 + wi 

02 = Al — 2A2 + V72, 

and hence the secular apsidal resonance variable 

#SAR = - W2 = 0i - 8 2 , 



(1) 

(2) 



(3) 



librate about 0°, which mean that the periapses are nearly aligned and that conjunctions of the 
planets occur when both planets are near periapse. In the above equations, Ai and A2 are the 
mean longitudes of the inner and outer planets, respectively, and Wj are the longitudes of peri- 
apse. There are three other systems wit 
Ferraz-Mello et al.ll2005l:lLee et al.ll2006l: 



73526 ( Tinnev et al 



2006; 



1 planets in 2:1 res onances: HD 82943 (IMayor et all 12004 ; 



Sandor et al. 



Beaugi_etal|l2008|), HD 128311 (|Vogt et al.ll2005h . and HD 
20071 ). although it should be noted that a pair of planets in 



1:1 re sonance is a plausible alterna tive for at leas t HD 82943 and H D 128311 (jGo zdzicwski & Konacki 



McArthur et al 



20061 ). In addition the HD 45364 dCorreia et~alll2009h . 55 Ca ncri ^M arcy et alJ booi: 
20041 ). HD 60532 (jPesort et alJbood : lLaskar fc Correialbood ). and HD 202206 |Cor"reia et al 



200 



5) 



systems have planets that are in 3:2, 3:1, 3:1 and 5:1 resonances, respectively. There are un- 
certainties in converting this data into the fraction of multiple-planet systems with mean-motion 
resonances. Some of the suspected resonant pairs may not be confirmed eventually (see, e.g., 
Fischer et al.ll2008l for 55 Cancri). On the other hand, the number of resonant pairs that remain 
undetected could be quite large, because the radial velocity variation due to two planets in reso- 
nance (in particular, 2:1) could be indistinguishable from that due to a single planet for certain 
planetary mas s ratio and orbital eccentrici t ies, given the precisio n levels of the existing radial ve- 
locity surveys (jAnglada-Escude et al.ll2008l : iGiuppone et al.ll2009i ). Nevertheless, the existing data 
indicate that ~ 20% of multiple-planet systems have mean- motion resonances. 

Mean-motion resonances can be easily established during planet formation by the convergent 
migration of planets due to interactions with the circumstellar gas disk. Two giant planets that 
are massive enough to open gaps in the disk individually can clear out the disk material between 
them rather quickly, and the outer planet is forced to migrate inward by the disk material out- 
side its orbit (and the inner p lanet outward if there is any disk material left inside its orbit) 
(IBryden et al.ll2000l : iKleyl l200(j ) . Both hydrodynamic and three-body simulations (with imposed 

y leads to capture 



migration for the latter) have shown that the convergence of t he orbits natural 



into mean-motion resonances (IBryden et al.l 120001: Kiev 



2002 



Nelson fc Papaloizoul 120021 : iPapaloizoul 120031 ; iThommes Lissauer 



2000 



Snellgrove et al 



2003 



2001 



Lee fc Peak 



Kiev et al.1 12004. 



2005; Lee 



20041 ). 
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The ubiquity of mean-motion resonances in extrasolar planetary systems and the ease of cap- 
ture into such resonances by convergent migrati on have prompted inve s tigations into the variety 



of stable mean-motion resonance configurations (ILee &: Peald [2002 



Ferraz-Mello et al 



2003 



Hadiidemetri ou fc Psycho 



oyc 

7E 



oJEo03; 



2003 



Ji et al 



Beauge et al 



2003 



Led 120041 ; IVovatzis & Hadiidemetrioul 120051 . 12009 ; iMarzari et al.l 120061 ; iMichtchenko et al 



2003 



2006 



Thommes fc Lissauer 



2006). The 2:1 mean-motion commensurability has received the most attention, because it is 
the most common one obs erved and includes the best case, GJ 876. With the exception of 
Thommes Sz Lissauer (j2003l ). all of the works just cited have focused on systems with two planets 
on coplanar orbits. For small orbital eccentricities, antisymmetric configurations with 9\ librating 
about 0° and 62 about 180° (as in the case of the Jovian satellites Io and Europa) are the only sta- 
ble 2:1 resonance configuration with both 9\ and #2 librating. For moderate to large eccentricities, 
the Io-Europa configuration is not stable, but there is a wide variety of other stable 2:1 resonance 
configurations, including symmetric configurations with both 6\ and #2 librating about 0° (as in the 
GJ 876 system), asymmetric configurations with 6\ and 82 librating about angles other than 0° and 
180° (some with intersecti ng orbits) , and antisymmetric configurations with 6\ ~ 180° and 62 ~ 0° 
(and intersecting orbits) 1*1 ILee (2004) has shown that the sequence of 2:1 resonance configurations 
that a system with initially coplanar and nearly circular orbits is driven through by continued mi- 
gration depends mainly on the planetary mass ratio m\/m2, if the migration rate is sufficiently slow. 
However, there are stable 2:1 resonance configurations (e.g., those with 9\ ~ 180° and 62 ~ 0°) 
that cannot be reached by the convergent migration of planets with constant masses and initially 
coplanar and nearly circular orbits. If real systems with these configurations are ever found, their 
origin would require a change in the planetary mass ratio m\/m2 during migration, multiple-planet 
scattering in crowded planetary systems, or a migration scenario involving inclination resonances 
QI2Qq3). 



Thommes fc Lissauerl (|2003l ) have studied the convergent migration of planets with non-coplanar 



orbits and found that, subsequent to the capture into the 2:1 eccentricity resonances, a capture into 
the 4:2 inclination resonances (which are the lowest order inclination resonances at the 2:1 com- 
mensurability) is possible if m\jm2 < 2. The 4:2 inclination-type mean-motion resonance variables 
are 



i>22 



2Ai 
2Ai 



4A 2 + 2^! 
4A 2 + 2Q 2 , 



(4) 
(5) 



where fij are the longitudes of the ascending node. The simultaneous librations of <j>\\ and $22 
mean that the mixed resonance variable 



Jia = 2Ai - 4A 2 + Oi + n 2 = {<f>n + </> 22 )/2 



(6) 



1 Throughout this paper, we often use "81 ~ x" as an abbreviation for "the libration of 8\_ about an angle x" (and 
similarly for the other resonance variables) when we describe a resonance configuration. 
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also librates. As the system enters the inclination resonances, the mutual inclination of the orbits 
can grow rapidly to tens of degrees. In some cases, the system eventually evolves out of the inclina- 
tion resonances, a nd the eccentricity resonance s switch to the Q\ ~ 180° and 62 ~ 0° configurations 
mentioned above. iThommes Lissauerl (|2003l ) showed an example with m\/m2 = 3, which does 
not have capture into the inclination resonances and remains nearly coplanar throughout its evo- 
lution, bu t we notice that the evolution of the eccentricities and 6j is different from that found by 
Led (|2004l ) when the outer orbital eccentricity e2 > 0.45. These configurations with > . 45 als o 
do not correspond to any of the other eccentricity resonance configurations found by iLed (|2004l ). 
As we shall see in £J3J they belong to a new family of 2:1 eccentrici ty re s onanc es that can be reached 
by migration if the migration rate is faster than that adopted by ILed (|2004l ) and m\/m-2 2. 



Gap-opening planets undergo type II migration on the disk viscous timescale, whose inverse is 

2 

(7) 



3v 
2^2 



9.4 x 10" 



n 



4 x 10 



-3 



H/a 



0.05 



(|Wardlll997l ). where a is the semimajor axis of the planet's orbit about a star of mass mo, a = da/dt, 
v = aH 2 Q, is the kinematic viscosity, a is the Shakura-Sunyaev viscosity parameter, H is the scale 
height of the disk, and P = 2-k/Q, ~ 27ra 3//2 /(Gmo) 1//2 is the orbital period. The uncertainties and 
radial variations in a and H/a mean that the migration rate can b e at l east a facto r of a few faster or 
slower than 10 _4 /P. Although both lThommes &: Lissauerl ( 20031 ) and ILed ( 2004! ) performed three- 
body simulations with imposed inward migration on the outer planet only, it is difficult to determine 
from the calculations in these papers that the different results at e2 > 0.45 for mj / 1712 = 3 are due t o 



differe nt migration rates, because they imposed migration in different ways. IThommes fc Lissauer 
( 2003 ) adopted H/a,2 oc a]/ 4 so that 02 is independent of 0,2 (with 0,2 = — 10 -5 AUyr^ 1 for most 
calculations) , and they imposed the migration in such a way that the migration does not slow 
down after the capture of an inner planet into resonance. ILed (120041 ) adopted constant H/a2 and 
performed calculations with 0,2/0,2 = — 10~ 6 /-P2 and — 10~ 4 /-P2> imposed in such a way that the 
migration slows down by a factor /3/{j3 + mi/m2), where (3 = 01/02 ~ 2 -2 / 3 , after the capture of 
an inner planet into 2:1 resonance. In this paper we examine systematically the effects of different 
migration rates (within the range expected for type II migration) using three-body integrations 
with migration imposed in the same way. 

We also examine systematically the effects of different eccentricity damping rates during mi- 
gration. Significant eccentricity damping can prevent the eccentricities from reaching high enough 
values for capture into the new 2:1 eccentricity resonances or the 4:2 inclination resonances. There 
is significant uncertainty in both the sign and the magnitude of the net effect of planet-disk in- 
teraction on the orbital eccentricity of the planet because of sensitivity to the distribution of disk 



material near the loca tions of the Lindblad and corotation resonances (jGoldreich fc Saril 120031 : 



Ogilvie fc Lubowl |2003| ). However, hydrodynamic simulations of two planets orbiting inside an 
outer disk have shown eccentricity damping of the outer planet, with K = ^2/^/^2/^ ~ 1 



(Kiev et al. 


2004, 


2005 


)• 



ber of simulations with eccentricity damping for non-coplanar and coplanar systems, respectively. 
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In particular, iThommes fc Lissauerl (j2003l ) have found that the critical value of K for capture into 
the 4:2 inclination resonances is between 2 and 5 for m\/m2 = 1. 

In £j2] we describe the numerical methods and initial conditions. In $3] we consider coplanar 
orbits and show that a resonant pair can evolve into a new fam ily of 2:1 eccentricity resonances 
if the migration rate is faster than that adopted by iLed ([2004) and m\jmi > 2, although the 
new family exists for m\/m,2 ^ 0.2. In §^] and [5] we consider non-coplanar orbits. We show that 
inclination excitation and capture into the 4:2 inclination resonances are p ossible for I mi 2 
(as we ll as m\/m2 < 2), if the migration rate is slower than that adopted bv lThommes Lissauer 
(|2003l ). and that the maximum value of K = 1 62 /e2 1 / 1 «2 / 0-2 1 for capture into the 4:2 inclination 
resonances is of the order of unity. Our conclusions are summarized and discussed in Sj6j 



2. NUMERICAL METHODS AND INITIAL CONDITIONS 

We consider systems consisting of a central star of mass m-o, an inner planet of mass mi, and an 
outer planet of mass m2, with mi/m.2 between 0.1 and 10. Unless stated otherwise, (mi+m2)/mo = 
10~ 3 . For the migration calculations starting with non-resonant orbits, the planets are initially on 
circular orbits, with the ratio of the orbital semimajor axes j3 = a\/o2 = 1/2 (far from the 2:1 
mean- motion commensurability where f3 ~ 2~ 2 / 3 ), and the outer planet is forced to migrate inward. 
The calculations presented in $3] are for configurations with exactly coplanar orbits, while those 
presented in §^1 and [5] are for configurations with initial mutual or bital inclination i mu = 0.01° , 



where the initial invariable plane is used as the z = reference plane. IThommes Sz Lissauerl ((2003) 
have found that the entry into the inclination resonances is not strongly influenced by the initial 
value of i mu , as long as it is < 1°. 

The t hree-body integrations with imposed migration are performed using the code described i n 



Lee) (J2004J), which is a modified version of the symplectic integrator SyMBA (IDuncan et al.l ll998). 
The outer planet is forced to migrate inward with a migration rate of the form 0,2/0,2 oc -P 2 _1 - 
The migration slows down by a factor /3/(/3 + mi/rri2), where (3 ~ 2~ 2 / 3 , after the capture of an 
inner planet into 2:1 resonance (see paragraph with eq. [7]). The input and output are in Jacobi 
orbital elements, and we apply the forced migration to the Jacobi 02 (and eccentricity damping 
to the Jacobi 62 for the calculations with eccentricity damping). To characterize the new family 
of 2:1 eccentricity resonances, there are also calculations in $3] with a change in mi / m? (and n o 



migration). The modified SyMBA code used for these calculations is also described in ILed ([2004]) 



3. A NEW FAMILY OF 2:1 ECCENTRICITY RESONANCES 



We begin with migration calculations of coplanar orbits without ecc entricity d amping. We 
consider mx/mi = 0.1, 0.3, 0.9, 1.0, 1.5, 2.3, 2.65, 3, 5, and 10 (same as in lLedl2004l ) and 0,2/02 = 
—0.5, —1, —2, —4, and —8 x 10~ 4 /i-2- Figures [T] and [2] show the evolution of the semimajor 
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axes dj, eccentricities ej, and eccentricity- type resonance variables Oj for the calculations with 
(mi + m-^/friQ = 10~ 3 , mi/mi = 3 and 02/02 = —0.5 x 10~ 4 /i-2 and —2 x 10~ 4 /P2) respectively. 
For 02/02 = —0.5 x 10~ 4 /P2 (Fig. [1]), the sequence of 2:1 resonance configurations after resonance 
capture — from (^1,^2) ~ (0°,180°) at small eccentricities to asymmetric librations of both 9\ 
and 6*2 at moderate to large eccentricities — is identical to that found by Lee (2004) for 1 02/02 1 < 
10~ 4 /i-2i but with larger libration amplitudes for faster migration rate. The system eventually 
becomes unstable at i/-P2,o = 3.96 x 10 4 , where i-2,0 is the initial outer orbital period. For the 
faster migration rate of 02/02 = —2 x 10~ 4 /i-2 (Fig. [2]), the sequence of resonance configurations 
is identical to that shown in Figure [T] (but with larger libration amplitudes) when e2 < 0.45 (and 
t/P2,o i$ 3000). However, the system enters a new family of 2:1 resonance configurations when 
e 2 ^ 0.45. The differences between the configurations in Figures [T] and [2] at e2 > 0.45 are most 
obvious in the plots of e\ and 9\. We confirm that the configurations in Figure [2] with e2 > 0.45 are 
stable resonance configurations by integrating the configurations at t/P^fi = 7000 and 10 4 forward 
with migration turned off and finding stable libration of 9j, with the centers and amplitudes of 
libration nearly identical to those just before the migration is turned off0 

Migration calculations with different mi/m,2 and 02/02 show that a system can enter the new 
family of 2:1 eccentricity resonances by fast migration if m\/m2 > 2. For (mi + 1712) /mo = 10 -3 , 
the transition occurs between 02/02 = — 1 x 10~ 4 /i-2 and —2 x 10~ 4 /i-2 for mi/m-2 = 2.65 and 3, 
and between 02/02 = —2 x 10~ 4 /i-2 and —4 x 10~ 4 /i-2 for = 2.3, 5, and 10. However, if 

02/02 is as fast as —8 x 10 _4 /i-2) the libration amplitudes are sufficiently large that the system 
becomes unstable soon after entering the new family. Calculations with twice the total planetary 
mass [(mi +m.2)/mo = 2 x 10~ 3 ] show that the critical migration rate for entry into the new family 
is roughly proportional to (mi + 7712) /mo- 

To find the small-libration-amplitude (or near exact resonance) counterpart for this new family 
and to determine the range of m\jm2 for which this family exists, we take the large-libration- 
amplitude configuration at t/P^jo = 7000 in Figure [2] and adjust the orbital parameters to obtain a 
small-libration-amplitude configuration with mi/m,2 = 3, e\ = 0.158, e2 = 0.702, Q\ = 1°, and 62 = 
98°. This small-libration-amplitude configuration is used as the starting point for two calculations 
in which m\jm2 is increased or decreased slowly to find a sequence of configurations with different 
mi/m2- The results are shown in FigureEJ with the calculations with dln(mi/m2) /dt = — 10 -6 /i-2,0 
and 10 _6 /P2,o along the positive and negative time axis, respectively. The inner eccentricity e\ 
increases (and outer eccentricity e2 decreases) with decreasing mi/m.2, and the system becomes 
unstable when mi/m-2 is decreased to about 0.2. The resonance configuration for a given mi/m,2 
from Figure [3] is then used as the starting point for slow inward (02/02 = — 10 _6 /P2) and outward 



Forced migration ca uses offsets in the libration centers of the resonance variables (|Leej 2004 



Murrav-Clav fc Chiang! [2005 ) . Both the offsets and the libration amplitudes increase with the migration rate. The 



offsets are typically much smaller than the libration amplitudes for the asymmetric configurations but could be no- 
ticeable for, e.g., the {61,62) ~ (0°, 180°) configuration that the system is first captured into (compare Figs. [T]and 
EJ). 
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(02/02 = 10~ 6 /P2) migration calculations to search for other resonance configurations with the 
same mx/m%. (The slow rate of change in mxlm-i or 02 in these calculations ensures that the 
libration amplitudes and offsets remain small.) Figure H]a shows the loci in the e\-ei plane of the 
stable resonance configurations from the migration calculations with m\jm% = 0.3, 1, 3, and 10. 
The initial conditions (dashed lines in Fig. [3]) are indicated by the triangles in Figure 0]ci, and the 
results from inward (outward) migration extend above (below) the triangles. In addition to the 
calculations shown in Figure El we perform several calculations in which different configurations 
along the locus shown in Figure H]a for mx/m^ = 0.3 are used as the starting point and m\jmi is 
decreased. In all cases, the system becomes unstable when mi/mj is decreased to about 0.2. Thus 
the new family of 2:1 eccentricity resonances does not appear to exist for m\/m2 < 0.2. 

The new family of 2:1 eccentricity resonances has e2 > 0.4-0.5, asymmet ric librations, and 
intersecting orbits, and it is distinct from any of the families found by iLed (J2004J). I n Figures \4$ )-d, 
we compare the new family (labeled IV) with the families (labeled I— III) found by [Leg (J200J) for 
mi/rri2 = 3, 1, and 0.3, respectively. Sequence I is the sequence reached by slow migration of planets 
with constant masses and initially nearly circular orbits; sequence II was found by a combination of 
calculations in which m\jmi is changed and slow migration calculations; and sequence III consists 
of configurations with (61,62) ~ (180°, 0°). For mi/m^ = 3 (Fig. 0]&), sequences I and IV come 
close to each other. In addition, like sequence IV, the confi guration s in sequence I with e2 > 0.34 
have intersecting orbits, as well as asymmetric librations (|Led 120041 ) . Thus it is possible to jump 
from sequence I to sequence IV if the libration amplitudes are large due to fast migration, as shown 
in Figure [2j For m\/m2 = 1 (Fig. (Us), sequences I and IV do not come close to each other. 
For m\/m2 = 0.3 (Fig. HH), although sequences I and IV come relatively close to each other at 
large e%, the configurations in sequence I at large e\ have (61,62) ~ (0°,0°) and non-intersecting 
orbits, and it is not possible to jump to sequence IV with asymmetric librations and intersecting 
orbits. Because the combination that sequences I and IV come close to each other and that some 
configurations along sequence I have asymmetric librations and intersecting orbits occurs only for 
m i/ m 2 ^ 2, we can understand why it is possible to enter the new family by fast migration only 
for mi/rri2 > 2, even though the new family exists for m\/m2 > 0.2. 



The existence of the new family of 2:1 eccentricity resonances was noted by lLee fc Thommes 



(|2004l ). IVoyatzis &: Hadjidemetrioul (|2005l ) have also discovered this family in their search for 
both stable and unstable asymmetric periodic orbits at the 2:1 resonance for three cases with 
m\/m2 = 0.54, 1, and 1.86. The stable periodic orbits are resonance configurations with zero 
libration amplit udes. In the range of m\/m,2 studied by IVovatzis fc Hadiidemetrioul (|2005l ) (i.e., 



m-\/m2 2.75. ILed I2004T) . s equen ce II has loci similar to those shown in Figures S]c and 0](i, and 



Vovatzis Hadiidemetrioul (|2005l ) found that sequences II and IV (the new family) are connected 
to each other by a sequence of unstable periodic orbits. 
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4. THE 4:2 INCLINATION RESONANCES 

We consider next migration calculations of non-coplanar orbits without eccentricity damping. 
We perform calculations with initial mutual orbital inclination i mu = 0.01° and migration rate 
^2/02 = —0.125, —0.25, —4 x 10 _4 /P2- Figures 03 [H and[7]show the results for m\lm% = 3 
and 02/122 = —2, —0.125, and —0.5 x 10~ 4 /P2 ; respectively. In each figure, we plot the evolution of 
the inclinations ij and the inclination- type resonance variables cftjj, as well as aj, ej, and 8j. The 
evolution of 0xi is nearly identical to that of (f>22, which means that f^i — 0,2 = (</>n — </>22)/2 (which 
is 180° initially due to our choice of the initial invariable plane as the z = reference plane) is close 
to 180° throughout and the ascending nodes are nearly antialigned. 

The fast migration calculation shown in Figure [5] is similar to that in Figure El but with non- 
coplanar orbits. The ej and 6j evolve as in the planar case shown in Figure [2] and enter the new 
family of eccentricity resonances when e2 > 0.45. There is no capture into the inclination resonances 
or excitation of the inclinations. Note, however, that the circulation of the inclination resonance 
variables changes from prograde to retrograde at about the same time as the entry into the new 



family of eccentricity resonances. The evolution in this figure is that found by lThommes &: Lissauer 



(2003) in their simulation with m\/m2 = 3. 



As in the planar case, the non-coplanar calculations with mi/mz = 3 and 0,2/0,2 slower than 
—2 x 10 _4 /P2 do not enter the new family of eccentricity resonances. However, unlike the planar 
case, the evolution for 0,2/02 slower than —0.5 x 10 _4 /P2 is qualitatively different from that for 
d 2 /a 2 = -0.5 and -1 x 10" 4 /P 2 - Fi gure [U shows a calculation with a slow migration rate (0,2 /02 = 
—0.125 x 10~ 4 /i-2)- The system is initially captured into the 2:1 eccentricity resonances only, and 
the initial evolution after capture is similar to that in the planar calculation with slow migration 
shown in Figure [TJ But starting at t/Pzp ~ 4.7 x 10 4 (when e2 ~ 0.45), the inclination resonance 
variables 4>jj change very slowly for about 6OOOP2.O) an d the inclinations increase rapidly. It is likely 
that this slow change of cftjj is associated with the proximity to the separatrix of the inclination 
resonances, since the circulation/libration period is infinite on the separatrix. We can understand 
qualitatively the almost exponential growth in the inclinations by noting that the lowest order 
inclination resonance terms at the 2:1 commensurability in the disturbing potential $ are second 
order and proportional to i\ cos <f>n, iiZ2Cos0i2, and i| cos 022- Thus the lowest order terms for 
dij/dt oc iJ 1 d&/dSlj are proportional to ijSixi(f)jj and ifcsin</>i2 (where k = 2 for j = 1 and vice 
versa), which can result in exponential growth if and $12 are not equal to 0° or 180° and 
change very slowly. At i/Pj.o ~ 5.3 x 10 4 , both fin and ^22 are captured into resonance and 
librate about 110°, and the inclinations increase slowly due to the continued migration forcing the 
system deeper into inclination resonances. We note that the simultaneous librations of 0j and <f>jj 
affect the values of 0j and ej during this phase (compare Figs. [Qand[6]). As in the case of the 
eccentricity resonances, asymmetric libration of 4>jj about an angle other than 0° or 180° is possible 
when the inclinations are not small and dij/dt is not dominated by the lowest order terms in the 
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disturbing potential]! We confirm that the configuration at, e.g., t/P^o = 8.0 x 10 4 in Figure [6] 
is indeed in stable inclination resonances by taking that configuration as the starting point for a 
three-body integration without forced migration and finding stable libration of (fijj about 110° and 
no secular change in ij throughout that integration. In Figure [6] the system eventually evolves out 
of the inclination resonances at t/P^Q ~ 1.0 x 10 5 , and t he eccentricity resonan c es sw itch to the 
9\ ~ 180° and 62 ~ 0° configuration. As mentioned in £pQ iThommes Sz Lissauerl (120031 ) have seen 
similar switching to the 6\ ~ 180° and 62 ~ 0° configuration in their simulations with mi/mj = 1. 

In contrast to the overall evolution on the migration timescale, the time spent in the phase 
with (fijj changing slowly and ij increasing rapidly is nearly independent of the migration rate for 
slow migration. Thus this phase takes up a larger and larger fraction of the total evolution time 
with increasing migration rate, and there is no longer a phase with (f)jj clearly in resonance if the 
migration rate 02/02 is as fast as —0.5 and —1 x 10~ 4 /P2- (Even faster migration rate would result 
in entry into the new family of eccentricity resonances and no inclination excitation, as discussed 
above.) For 0,2/(12 = —0.5 x 10~ 4 /i-2 (Fig. [7J), the rapid inclination excitation phase occurs from 
t/P2 ~ 1-2 x 10 4 to 2.4 x 10 4 . Then <pjj, as well as 0\, alternate between libration and circulation 
for about 6000P2,o> before (f>jj change to circulation only and the eccentricity resonances to the 
9\ ~ 180° and 82 ~ 0° configuration, with 6\ nearly circulating but spending most of its time 
around 180°. The oscillations of the inclination resonance variables cj>jj between t/Pifs ~ 1-2 x 10 4 
to 2.4 x 10 4 in Figure [7J might lead one to think that <fijj are in resonance and librating about 
equilibrium values and that the rapid increase in the inclinations is due to continued migration 
forcing the system deeper into inclination resonances. However, this would be inconsistent with 
our earlier observation that the duration of this phase is nearly independent of the migration rate 
for slow migration. To show that this rapid inclination excitation is in fact not due to migration 
forcing, we take the configuration at t/i"2,o = 2.0 x 10 4 in Figure [7J as the starting point for a 
three-body integration without forced migration. The results are shown in Figure [8] As we can 
see, the inclinations continue to increase rapidly for about 4000i-2,o even without forced migration. 
Furthermore, the evolution of all the plotted variables for the first 10 4 P2,o i n Figure [8] without 
migration is similar to that between t/P^fl = 2.0 x 10 4 and 3.0 x 10 4 in Figure [7J with migration. 
Figure [8] also shows us what would happen if the migration stops due to, e.g., disk dispersal when 
the system is in the phase with rapid inclination excitation. The inclinations would continue to 
increase for a while, and the inclination resonance variables would eventually end up in large- 
amplitude libration (alternating with circulation to varying degree). 

Figure [9] summarizes the results for (mi + 7712) /mo = 10 -3 and different ra\jvri2 and 0,2/02- For 



3 In the limit of the circular, planar, restricted, three-body problem, one can identify the terms that give rise 
to asymmetric libration for the n:l (not just 2:1) exterior resonance as coming from the indirect part of the dis- 
turbing potential, and there is a qu alitative physical explanation based on the indir ect acceleration imparted on the 
test particle over a synodic period jPan fc Saril 120041 : IMurrav-Clav fc Chianell2005l) . This type of analysis has not 
been generalized to either the planar two-planet problem with two resonance variables 6\ and 62 or the inclination 



resonances. 
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m\/m2 > 2 and fast migration (the region labeled E in Fig. [9]), the eccentricity resonances enter 
the new family, and there is no capture into inclination resonances or excitation of the inclinations. 
For slow migration (the region below the solid line in Fig. [9]), the inclination resonance variables 
are captured into libration after a phase with fyj changing slowly and ij increasing rapidly. 
The inclination resonance configur ation is symmetric with 4>jj librating about 180° in the region 



labeled S with m\/rri2 % 2.5 (see iThommes &: Lissauerl 120031 for an example with m,\jm% = 1), 
and it is asymmetric with <f>jj librating about an angle other than 0° or 180° in the region labeled 
A with mi/ m2 > 2.5 (e.g., Fig. [6]). For intermediate migration rate (and also fast migration 
rate if mi/m-i < 2), we typically see the rapid inclination excitation phase, but not a phase with 
4>jj clearly in resonance (e.g., Fig. [7J). The inclination excitation can be partial, with the mutual 
inclination reaching a maximum of ~ 1° or less, if mi/m2 is large (in particular mi/m2 = 10). 

In the phase with simultaneous librations of the eccentricity and inclination resonance variables, 
we can see from the definitions of 0j (eq. [TJ [2]) and fyj (eq. [4]-[5]) that the arguments of periapse 
ujj = Wj — = 9j — <pjj/2 also librate. For mi/m2 < 2.5 with symmetric libration of <pjj, tOj librate 
about ±90° (i.e., the periapse is on average 90° ahead of or behind the ascending node), while for 
mi/m2 > 2.5 with asymmetric libration of 4>jj, the libration of Uj is also asymmetric. 



5. EFFECTS OF ECCENTRICITY DAMPING 

As we mentioned in fJJJ sufficient eccentricity damping can prevent the eccentricities from reach- 
ing high enough values for inclination excitation and/or capture into the inclination resonances. 
In order to study the effects of eccentricity damping, we repeat the non-coplanar calculations in 
$4] with the ratio of eccentricity damping to migration of the outer planet, K = \ez/ei\/\(i2/ az\, 
ranging from 0.25 to 8. 

We consider first the calculations with slow migration (0,2/0,2 below the solid line in Fig. 
[9]). Figure [10] shows the evolution of the mutual inclination i mu for mi/m,2 = 0.3, 1.5, and 5.0, 
02/02 = —0.125 x 10 _4 /P2i an d different K. As K increases from zero, the system enters the rapid 
inclination excitation phase and the subsequent capture into the 4:2 inclination resonances later 
and later, because the eccentricities grow slower and slower. However, when K exceeds a critical 
value, the eccentricities never reach high enough values for inclination excitation and capture into 
inclination resonances. The critical value of K is ~ 1.4 for mi/m2 < 0.3, ~ 2.8 for mi/7712 ~ 0.9— 
1.5, and « 0.7 for m 1 /m 2 > 2.65 (Fig. ITT]) . 

For faster migration rate, the effects of eccentricity damping on the evolution of the system 
can be more complicated. For example, the calculation shown in Figure [12] is similar to that 
in Figure [7] (mi/m2 = 3 and 0,2/02 = —0.5 x 10 _4 /P2) but with K = 0.25. In this case, the 
eccentricity damping results in clear libration of the inclination resonance variables <pjj after the 
rapid inclination excitation phase. Nevertheless, the critical value of K as a function of mi/m2 
shown in Figure [TT] also summarizes the results for migration rate up to d 2 /a 2 = -2 x 10" 4 /P 2 , if 
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lt is interpreted as the critical value for inclination excitation, which may or may not be followed 
by a phase with (fijj clearly in resonance. For CI2/02 = —4 x 10 -4 /P2 (the maximum migration rate 
studied), the critical value of K is modified at large mi/m.2, with none of the calculations with 
vn\jm2 > 5 showing inclination excitation. 



6. CONCLUSIONS 

We have investigated the effects of different migration rates on the capture into and evolution in 
eccentricity and inclination resonances at the 2:1 mean- motion commensurability by the convergent 
migration of two planets. We focused on systems with orbits that are initially slightly inclined with 
respect t o each othe r. The system is first captured into the sequence I of 2:1 eccentricity resonances 



found bv lLed (|2004l ). the same as in the case of exactly coplanar orbits. If the migration rate is fast 
and m\lm% > 2, the subsequent evolution is also identical to the coplanar case, with the eccentricity 
resonances entering a new family (sequence IV), and there is no inclination excitation or capture 
into inclination resonances. The new family of 2:1 eccentricity resonances (with e2 > 0.4-0.5, 
asymmetric librations, and orbits that intersect if they are exactly coplanar) exists for mxjm^ > 0.2, 
but it is possible to evolve into this family by fast migration only for m\jm2 > 2. If the migration 
rate is slow, the system subsequently enters a phase with the 4:2 inclination resonance variables 
(j)jj changing slowly and the inclinations increasing rapidly, before it is captured into 4:2 inclination 
resonances. The inclination resonance configuration is symmetric, with 0xi ~ 022 ~ 180°, if 
m -i/ m -2 ^ 2.5 and asymmetric if m\jra2 > 2.5. For intermediate migration rate (and fast migration 
rate if m\jm2 < 2), there is typically a rapid inclination excitation phase, but not a phase with (ftjj 
clearly in resonance. We have also studied the effects of different eccentricity damping rates during 
migration and found that the maximum value of K = \&2l &2\l\&2l &2\ for inclination excitation 
(which may or may not be followed by a phase with cfijj clearly in resonance if the migration rate 
is not slow) ranges from ~ 0.7 for m\/m2 > 2.65 to ~ 2.8 for mi/771,2 ~ 1. Since the evolution is 
sensitive to the rates of migration and eccentricity damping within the ranges expected for type II 
migration due to planet-disk interaction, the discovery of extrasolar planetary systems with certain 
combinations of mass ratio and 2:1 resonance geometry would place a constraint on the strength 
of eccentricity damping during migration, as well as on the rate of migration itself. 

There are several effects of disk-planet interaction that were neglected in our analysis and may 
require further investigations. We have focused on inward migration and eccen t ricity damping o f 



the outer planet, because previous hydrodynamic simulations (e.g., iKleyl 120001 : iKley et al.1 12004 ) 



have shown that the disk inside the inner planet's orbit , not just the disk material between the 



planets, should be cleared rapidly. However, ICrida et al.1 (120081 ) have recently shown that a better 



numerical treatment of the inner disk may result in a slower depletion of the inner disk and that the 
eccentricity damping from the inner disk could be important in explaining the observed eccentricities 
of resonant pairs such as that in the GJ 876 system. On the other hand, when the nearby disk 
mass is comparable to the planet mass (i.e., in older, partially depleted disks), planets will undergo 
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type II migration at significantly less than the disk's viscous advection speed (jSver &: Clarkdll995l ). 
so that the inner disk "outruns" the planet and eve ntually leaves an inner hole, no matter how 
small the inner boundary radius. The simulations of iThommes et ajj (|2008l ) suggest the majority 
of planets form late enough in their parent disk's lifetime that such holes are ubiquitous. 

Planet-disk interaction can also affect the orbital inclination of a planet. The net effect of 
inclination damping by secular interactions and excitation by interactions at mean-motion reso- 
nances depends on the disk parameters, but any net damping should be o n a timescale comparable 
to or longer than the migration timescale (e.g., iLubow Qgilvid l200ll ). This is likely too slow 
to affect the rapid inclination excitation phase, but may result in equilibrium inclinations if the 
system is subsequently captured into inclination resonances and the inclinations are excited slowly 
by continued migration. 

We have also neglected the secular apsidal and nodal precessions induced by the disk, which 
could change the sequence of resonance capture by changing and splitting the locations of the 
various resonances at the same mean-motion commensurability. Kiev et al.l ((2005) have performed 
coplanar three-body integrations of the GJ 876 resonant pair with additional apsidal precession 
and found that the eccentricity resonances 9\ and 62 are captured into libration in a sequence 
that differs little in order or timing from the case without additional apsidal precession. This 
can be explained by the fact that the 2:1 eccentricity resonances are first order, which means 
that the resonance-induced retrograde apsidal precession is proportional to 1/ej and much larger in 
magnitude than the disk-induced prograde precession for small e,-. On the other hand, disk-induced 
nodal precession could have a larger effect, because the 4:2 inclination reson ances are second order 
and th e resonance-induced nodal precession is roughly constant for small ij. IThommes fc Lissauer 
(|2003l ) have performed some non-coplanar calculations with additional apsidal and nodal precessions 
and did not find any significant difference from the calculations without additional precessions. 
Although the adopted disk surface density is 5 times that of the minimum mass solar nebula, they 
assumed an outer disk with an inner edge that is likely too far (20 Hill radii) from the outer planet's 
orbit, and the amount of precession induced by the disk is determined primarily by the material 
closest to the planet. 



Adams et al.l (|2008l ) and lLecoanet et al.l (|2009l ) have recently examined the effects of turbulence 



in circumstellar gas disks on mean-motion resonances in extrasolar planetary systems. They have 
found that stochastic perturbations due to turbulence could prevent planets from staying in resonant 
configurations and that planetary systems with mean-motion resonances should be rare. This 
appears to be inconsistent with the observational evidence discussed in £[TJ One possible explanation 
is that these studies assumed full magnetorotational turbulence, whereas circumstellar disk models 
usually exhibit an extensive dead zone around the midplane, where the ionization fraction is low 



and the disk is magnetorotationally stable due to ohmic dissipation (e.g.. lGammidll996l : ISano et al 



2000: lTurner et alJl2007l : lllgner &: Nelsonll2008l ). For weak turbulence, the turbulence may generate 
larger libration amplitudes than in smooth migration and allow, e.g., the jump from sequence I to 
sequence IV to occur at a slower migration rate. 
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Finally, a better understanding of the capture into the inclination resonances is needed. Since 
the time spent in the phase with (j)jj changing slowly and ij increasing rapidly is nearly independent 
of the migration rate for slow migration (see in the limit of very slow migration, there is an 
almost instantaneous jump in the inclinations at the time of the inclination resonance capture, if 
we measure time in units of the migration timescale. This is not what one would expect if the 
capture into the inclination resonances can be modeled by the appearance of additional equilib- 
rium j3oints_arid_sej^^^ the Hamiltonian theory of a single second-order resonance (see, 



e.g., Murray Sz Dermottlll999l ). The Hamiltonian approach is based on the assumption that each 
resonance is encountered individually, which is clearly not the case in our problem. In particular, 
the system is already in eccentricity resonances when it enters the inclination resonances. If the 
system is captured into just an inclination resonance, the capture requires Clj « — Ai + 2A2. On the 
other hand, if the system is already in eccentricity resonances so that Wj ~ — Ai + 2A2, then the 
capture into the inclination resonances requires ujj ~ 0. 

It is a pleasure to thank Stan Peale for informative discussions. This research was supported 
in part by NASA grant NNG06GF42G (M.H.L.) and a grant from NSERC Canada (E.W.T.). 
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t/P 2 ,o 

Fig. 1. — Evolution of the semimajor axes di and 02, eccentricities e\ and C2, and 2:1 eccentricity-type 
mean-motion resonance variables 6\ = Ai — 2A2 + roi and 02 = Ai — 2A2 + VJ2 for a differential migration 
calculation of coplanar orbits without eccentricity damping. The mass ratios (mi + •m-^jmo = 10~ 3 and 
mi/rri2 — 3. The outer planet is forced to migrate inward with 0,2/0,2 — —0.5 x 10~ 4 /P2- The semimajor 
axes and time are in units of the initial orbital semimajor axis, 02, 0, and period, i-2,0 °f the outer planet, 
respectively. The sequence of resonance configurations after resonance capture — (^1,6*2) ~ (0°,180°) — > 
asymmetric librations — is identical for |d 2 /a 2 | < 10~ 4 /P 2 , but with larger libration amplitudes for faster 
migration rate. 




t/P M 

Fig. 2. — Same as Fig. [1] but for the faster migration rate of 02/02 = —2 x 10~ 4 /P2- The system enters a 
new family of 2:1 resonance configurations when e2 > 0.45 (and t/P2 } o > 3000). 
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Fig. 3. — Evolution of the eccentricities ei and e2, eccentricity- type resonance variables #i and 02, and 
mass ratio mi/ni2 for calculations in which a configuration in the new family of 2:1 eccentricity reso- 
nances with m\/m2 = 3 is used as the starting point and m\jm,2 is increased and decreased. The start- 
ing configuration with small libration amplitudes is obtained by adjusting the orbital parameters of the 
large-libration-amplitude configuration at t/Pz^ = 7000 in Fig. [2] The results from the calculations with 
d\n(nii/m2)/dt = — 10~ 6 /P2,o and 10 _6 /P2,o are plotted along the positive and negative time axis, re- 
spectively. The system becomes unstable when mi/m 2 is decreased to about 0.2. The configurations with 
myjmi = 0.3, 1, 3, and 10, indicated by the dashed lines, are used as initial conditions for calculations in 
Fig. Hi. 
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Fig. 4. — Loci in the ei-e% plane of coplanar 2:1 resonance configurations, (a) Configurations in the new 
family of eccentricity resonances from inward and outward migration calculations with initial conditions 
{triangles, oriented to indicate the direction for inward migration) from Fig. [3] for mx/m<2 — 0.3 , 1, 3, and 
10. Comparison of the new family (labeled IV) with the families (labeled I— III) found by I Lee! ( 2004 ) for 
m\/m2 = (b) 3, (c) 1, and (d) 0.3, respectively. It is possible to enter the new family by fast migration only 
for nii/ni2 > 2 because the combination that sequences I and IV come close to each other and that some 
configurations along sequence I have asymmetric librations and intersecting orbits occurs only for mi/m 2 > 2 
(see text for details). 
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Fig. 5. — Evolution of the semimajor axes a\ and 02, eccentricities ei and e2, 2:1 eccentricity-type resonance 
variables #1 and 02, inclinations i\ and 12, and 4:2 inclination-type resonance variables </>n = 2Ai — 4A2 + 2f2i 
and ^22 = 2Ai — 4A2 + 2f^2 for a differential migration calculation of non-coplanar orbits without eccentricity 
damping. The mass ratios (mi +TO2)/too = 10 -3 and m\jm2 — 3, and the initial mutual orbital inclination 
«mu = 0.01°. The outer planet is forced to migrate inward with the fast migration rate of 0,2/0,2 — —2 x 
1CP 4 /P2- The eccentricity resonances enter the new family as in the planar case shown in Fig. [2J and there 
is no capture into the inclination resonances or excitation of the inclinations. 
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Fig. 6. — Same as Fig. [5j but for the slow migration rate of 0.2/02 = —0.125 x 1CP 4 /P2- The system 
is initially captured into the 2:1 eccentricity resonances only. There is a phase from t/P^fi ~ 4.7 x 10 4 to 
5.3 x 10 4 with <j>jj changing slowly and ij increasing rapidly before <pjj are captured into libration. The system 
eventually evolves out of the inclination resonances at i/i-^.o ~ 1-0 X 10 5 , and the eccentricity resonances 
switch to the 6\ « 180° and #2 ~ 0° configuration. 
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Fig. 7. — Same as Fig. El but for the intermediate migration rate of 02/02 = —0.5 x 10~ 4 /P2- Fhe rapid 
inclination excitation phase occurs from t/P%fi « 1.2 x 10 4 to 2.4 x 10 4 . Then <pjj and 0\ alternate between 
libration and circulation for about 6OOOP2.0, before (f>jj change to circulation only and the eccentricity 
resonances to the 6>i « 180° and 9 2 « 0° configuration. 
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Fig. 8. — Evolution for the calculation in which the configuration at t/P2,o = 2.0 x 10 4 in Fig. [7] is used 
as the starting point for a three-body integration without forced migration. The inclinations continue to 
increase rapidly for about 4000P2,o, and the evolution of all the plotted variables for the first 10 4 P2,o is 
similar to that between t/P 2 ,o = 2 -0 X 10 4 and 3.0 x 10 4 in Fig. [7]with migration. 
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Fig. 9. — Types of evolution for different mi/m 2 and a 2 /a 2 . The results are from migration calculations 
with (mi + JTi2)/mo = 10~ 3 , non-coplanar orbits, and no eccentricity damping. In the region labeled E, the 
eccentricity resonances enter the new family, and there is no capture into inclination resonances or excitation 
of the inclinations. In the region below the solid line, the inclination resonance variables <fijj are captured 
into libration (symmetric in the region labeled S and asymmetric in the region labeled A) after a phase with 
(fijj changing slowly and ij increasing rapidly. In the unlabeled region, there is typically a rapid inclination 
excitation phase, but not a phase with <pjj clearly in resonance. 
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Fig. 10. — Evolution of the mutual inclination i mu for toi/to 2 = 0.3, 1.5, and 5.0, a 2 /a 2 = — 0.125 x 1CP 4 /P 2 , 
and different eccentricity damping ratio K = |e 2 /e 2 |/|a 2 /a 2 |. The eccentricities never reach high enough 
values for inclination excitation and capture into inclination resonances when K exceeds a critical value. 
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Fig. 11. — Critical value of K as a function of mi/m^. The critical value is for capture into inclination 
resonances for slow migration (02/02 below the solid line in Fig. [9]) and for inclination excitation (which 
may or may not be followed by a phase with cj)jj clearly in resonance) for migration rate up to 02/02 = 
-2 x 10~ 4 /P 2 - 




Fig. 12. — Same as Fig. [7j but with eccentricity damping ratio K = 0.25. In this case, the eccentricity 
damping results in clear libration of 4>jj after the rapid inclination excitation phase. 



